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Abstract 

The superfluid properties of a two-state Fermi mixture in an optical lattice are profoundly 
modified when an imbalance in the population of the two states is present. We present analytical 
solutions for the free energy, and for the gap and number equations in the saddle-point approxima- 
tion describing resonant superfluidity in the quasi-two-dimensional gas. Inhomogeneities due to the 
trapping potentials can be taken into account using the local density approximation. Analyzing the 
free energy in this approximation, we find that phase separation occurs in the layers. The phase 
diagram of the superfluid and normal phases is derived and analytical expressions for the phase 
lines are presented. We complete the investigation by accounting for effects beyond mean-field in 
the BEC limit where the system is more properly described as a Bose-Fermi mixture of atoms and 
molecules. 



The formation of pairing correlations in a mixture of two types of fermions is frustrated 
when the number of fermions in each state is unequal. The effect of such "spin imbalance" on 
pairing has been investigated theoretically since Clogston's seminal work [lj|for conventional 
superconductors, and has led to the prediction of neve, pairing states fl, ^ -levant 
for color superconductivity in dense quark matter |3| and for neutron-proton pairing in 
asymmetric nuclear matter J]. The experimental study of Fermi superffuids with imbalanced 
spin population has on.y recently become possibfe in u.traco.d atomic gases flQ, renewing 
theoretical interest in imbalanced superfluidity [3, [8|, [9|. 

In these experiments, the interaction strength between ultracold fermions can be precisely 
controlled through the use of Feshbach resonances. These scattering resonances allow to tune 
the s-wave scattering length from a large negative value, giving rise to a Bardeen-Cooper- 
Schrieffer (BCS) superffuid [10|], to a large positive value where a Bose-Einstein condensate 
(BEC) of weakly bound molecules is formed [ll]. 

Also the trapping geometry confining the Fermi gas can be precisely controlled experi- 
mentally. Of particular interest is the possibility to impose a crystalline potential through 
the use of optical lattices. These allow to experimentally mimic theoretical lattice models 
such as the Hubbard model. One-dimensional optical lattices allow to study stacked layers of 
superfluid and create a geometry analogous to layered (cuprate) high-temperature supercon- 
ductors. Optical lattices have been used to demonstrate superfluid behavior of condensates 



[12], and to probe the Mott-superfluid transition 131 ] . So far, fermionic superfluidity in an 
optical lattice has only been studied [yj] with balanced Fermi gases. 

We investigate the effect of spin imbalance on the superfluid properties of a Fermi gas in 
an optical lattice. In particular, we will examine the case of a one- dimensional optical lattice 
generated by two counterpropagating laser beams (parallel to the z-axis) with wave length A. 
These laser beams generate a periodic potential V sin 2 {2nz/\) .When loaded in this optical 
lattice, the gas forms a stack of typically a few hundred quasi-2D layers containing several 
thousands of atoms each. The interaction between atoms within a given quasi-2D layer can 
be modelled by a 2D contact interaction whose strength g depends on the model cutoff K c 



and on the energy of the scattering atoms through 

In (E/E b ] 
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U (2ir)3 («)Vm-B + fe' 
Here, m is the mass of the atoms and Ej, is the energy of the bound state that always exists 



2 



in two dimensions, given by 

£t = ^exp(WiY (2) 
7T \ a s J 

with a s the (3D) s-wave scattering length of the fermionic atoms, ujl = a/87t 2 Vo/ (mA 2 ) and 
£ L = ^h/{muo L ) and C « 0.915 (cf. Ref. Q). 

In this letter we derive an analytical expression for the free energy of the gas with density 
n in a layer in the optical potential. Extremizing this free energy with respect to the 
superfluid gap allows to set up and solve the gap equation. We derive analytical expression 
for both the gap and the chemical potential, for fixed imbalance. However, we will argue 
that the imbalanced gas may be unstable with respect to phase separation into a balanced 

fj 

superfluid and a halo of excess carriers, similar to the three-dimensional case [8j], and we 
derive the phase diagram for the gas, again retrieving analytic expressions for the phase 
boundaries. 

The partition sum of the quasi-2D Fermi gas can be written as a path integral over the 
exponential of the action functional S for the fermionic fields ipk,<T,ipk,<T where a denotes 
the spin. We write k = {k, u n } for the 2D wave number k and the Matsubara frequency 
u n = (2n + l)vr//5 where (3 = 1/(/cbT) is the inverse temperature: 

Z = J VipVtPexp{-S}, (3) 

with 

S = i ^2 [ i7T^ \ E </V [-w n + k 2 - /v] 4>k,cr + ffipk,^-k,iip-k,iipkA \ ■ ( 4 ) 

Here, we use units such that h = kp = Ep = 1. The number of spin-up and spin-down 
fermions is controlled through the chemical potentials /x-f and The partition sum cor- 
responding to the action functional (j4]) can be calculated following the standard procedure 
of introducing the Hubbard- Stratonovic decomposition, integrating out the Grassman vari- 



ables, and performing the saddle-point approximation |17l . 1 1 8f | with a constant gap A. We 
therefore neglect the possibility of two-dimensional FFLO type states with a modulated 
order parameter 0, It is useful to express the results as a function of \i = (/i-f + //j.)/2 
determining the total number of fermions and £ = (//-[■ — jt/j_)/2 expressing the imbalance in 
chemical potentials. We find Z sp = exp{— j3VL sp ] with 

Q sp = -— - \ f {log [2 cosh(/3C) + 2 cosh((3E k )] - (3£ k } . (5) 

g [3 J (2ny 



Here the fc-integration is convergent at large wave vectors so that the wave vector cutoff 
K c can be sent to infinity, A represents the pairing gap, = k 2 — fx is the free particle 
spectrum, and = a/£| + A 2 is the Bogoliubov spectrum. The saddle-point free energy 
Q sp allows to derive the gap and number equations. A similar result was derived for the 
homogeneous 3D gas by Iskin and Sa de Melo [3]. For the 2D case however, we succeeded 
to perform the wave-number integrations analytically in the limit of temperature going to 
zero: 

A 2 



8nn sp = -A 2 In (E b ) - fx (vV + A 2 + //) + A 2 In (a/V + A 2 - /i) - — 
- 2C (H - kl) + + A 2 - £a V 7 !^ + A 2 In (^=^| 

where £ a fi = k 2 b — fx and for £ > a/ A 2 + (min[/i, 0]) 2 

fc^^ /max^-y/^ -A^O] 



(6) 
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For £ a/A 2 + (min[/i, 0]) 2 the terms containing k a and k b in ([6]) vanish. The Bogoliubov 
energy E^ for quasiparticle excitations with wave vector k = k a , k b becomes equal to the 
imbalance of the chemical potentials (. All quasiparticle states with k e [k a , kb] have a 
Bogoliubov energy less than (, and contribute as normal state particles, rather than as 
superfluid excitations, to the free energy. 

Iskin and Sa de Melo j^] solve the gap and number equations for the imbalanced three- 
dimensional superfluid numerically with respect to (,fx and A for a given imbalance 5n/n 
and a given interaction strength g. In the two-dimensional case, this procedure applied to 
(EJ), yields analytical results, summarized in Table [B However, the question can be raised 
as to whether the state with fixed 6n/n is stable with respect to phase separation into a 
phase with a balanced superfluid (Sn/n = 0) in the center of the trap and an imbalanced 
halo of excess spin component surrounding it. Indeed, in the experiment of Ketterle and 



co-workers [20(] for a three-dimensional gas with a nonuniform trapping potential, the excess 
spin component is expelled in a shell surrounding a spin-balanced superfluid . Will this be 
similar for a 2D gas within a layer of the optical potential ? 

If we consider a 2D trapping potential V (r) that varies slowly on the length scales set by 
the interparticle distance £ = n" 1 ^ 2 and the superfluid gap = a/ h 2 /mA, we can apply the 



A 2 



E b /2 < (5n/nf E F (Sn/n) 2 E F < E b /2 
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2E b E F [1-h (E F /E b ) — 

\ n 

E F E F -^(l-h (E F /E b ) *£\ 



TABLE I: Analytical solutions for A 2 and \i at fixed imbalance 5n/n and binding energy E b , where 
h(x) = max [\/2x, l] . For E b /(2E F ) < (5n/n) 2 superfluidity is suppressed. There is a qualitative 
change in the dependence of A 2 and \i on E b when E b becomes equal to 2E F ; this value can be 
intepreted as separating a 'weak coupling' from a 'strong coupling' regime. 



procedure outlined by De Silva and Mueller 21( for the imbalanced gas in a regular three 



dimensional trap. That is, we calculate n (r) and Sn (r) in the local density approximation, 
where the local average chemical potential equals fi = /io — V (r) and the difference between 
the chemical potentials C = Co is constant in space for a spin-independent trapping potential: 

n (r) = - dn sp /d^\^_ vm=<o (8) 
5n(r) = -dn sp /dC\^ _ v{rU=(0 . (9) 

The gap is then found by extremizing Q sp with respect to A, for /i and ( fixed by the local 
density approximation, rather than for 5n/n fixed. The validity of this procedure depends 



on whether the loca 
experimental setups 



density approximation is justified; this need not be the case for all 

3. 



For ( = 0, we find that, the free energy Q sp shows a single minimum at Abai = 
\/2E b (/i + E b /2) [l^], where its value is fibai = — + E b /2) 2 /(4tt). This minimum rep- 
resents the superfluid state (since A ^ 0). As the imbalance ( is increased, Q sp develops 
a second minimum around A = 0, representing the normal state. The free energy of the 
normal state is fi = — [(/i + C) 2 ©(Ai + C > 0) - C) 2 B(/x - C > 0)] /(8tt). We have verified 
numerically that no other minima of Q sp except the ones at A = Ab a i and A = occur. 

Upon increasing ( the free energy curves are only affected in the region A < (. The 
minimum at A = Abai, representing the superfluid state, is therefore not affected by changes 
in ( until it ceases to exist: flh&i is independent of (. From this it follows that the superfluid 
state does not sustain imbalance (5n = since dQt, a \/d( = 0). The free energy of the normal 
state, fl , does depend on (, so that the normal state supports imbalance, as expected. This 
scenario is similar to that for a 3D imbalanced Fermi mixture 
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FIG. 1: The phase diagram for the imbalanced (two-dimensional) Fermi gas is shown as a function 
of the average chemical potential n and the difference between the chemical potentials Three 
phases can be identified: a balanced superfluid (SF), an imbalanced normal sate (N) and a fully 
polarized normal state (NP). The arrow indicates the path in the phase diagram that is traversed 
when moving away from the center of a 2D trap towards the edges. The corresponding overall 
shell structure is illustrated in the inset. At very strong coupling (large E b ) we expect corrections 
beyond mean field to push the phase line down into the shaded region, where also a mixture of 
bosonic molecules and fermionic atoms can appear. 

To evaluate the phase boundary between the balanced superfluid and the imbalanced 
normal state, we should make a distinction between the cases fi > ( and fi < (. For 
fi > C the thermodynamic potentials of the superfluid and normal state are equal for fj,Eb = 
C 2 — (E b /2) 2 , whereas for ( > fi > they are equal for y/2E b )/(\/2 — 1). In the 

normal state, the chemical potential of at least one of the spin components (up or down) 
should be positive for there to be any particles at all. This means that /i + C > or /x — ( > 0. 
For C > 0, the first condition is easier to fulfill, so that it determines the boundary in the 
phase diagram with the "empty phase" (no particles). 

Figure [1] shows the phase diagram of the superfluid (SF), normal (N), fully polarized 
normal (NP) and empty phases as a function of [i/E b and (/E b . This phase diagram looks 
qualitatively similar to the one obtained in 3D by mean field theory in Ref. j^ll ]. 




Within this local density approximation, moving from the center of the trap towards the 
edge corresponds to moving down along the arrow in figure [H From the center to the edge 
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we encounter first the balanced superfluid (SF), then a shell of spin-imbalanced normal gas, 
and finally a shell of fully polarized normal gas, as illustrated in the inset. The superfluid 
state in this treatment is balanced in the limit of temperature zero, but allows for imbalance 
for any T > 0. 

In a typical experiment with a ID optical lattice, several hundreds of layers are present. 
As long as the Fermi energy and the superfluid gap are much larger than the tunneling rate 
between the layers, the pairing is essentially a phenomenon that takes place within each 
layer individually [22J and each layer can be treated in the local density approximation by 
setting fiQ = fi — where is the trapping potential at the center of the ith layer. 
The tunneling between the different layers is between layers of zero imbalance, so that the 



results of Refs. 



13, 



23| | can be used. 

The preceding saddle point calculation implicitly makes the mean field assumption that 
the typical size of the paired state l& is much larger than the distance between the fermionic 
atoms £. Deep in the BEC limit, the system is in fact more suitably described as a Bose- 
Fermi mixture of strongly bound bosonic pairs and excess majority component atoms. We 
again can treat the 2D Bose-Fermi mixtures in the fi-( plane analytically. The energy density 
of the Bose-Fermi mixture is given by E = 2im 2 p + Qbb^bI^ + Qbf^b^f where the density 
of bosonic molecules n# and of fermionic atoms rip are positive. The coupling constant gss 



24|, but not that for 



for dimer-dimer scattering in 2D has been calculated by Petrov et al. 
dimer-atom scattering, gsF- However in the limit of very low energy scattering, dimensional 
arguments show that to first order the dimer-dimer and dimer-atom scattering amplitudes 
are equal. The chemical potentials are fip — Airnp + Qbf^b for the fermionic atoms and 
A*B — 9bb u b + Qbf^f for the bosonic molecules. 

A minimal energy state has a positive curvature as a function of the densities, i.e. the 
eigenvalues of = d 2 E / dnidrij have to be positive. A necessary condition is that its 
determinant is positive: ^rqbb — 9bf > 0- Inverting the expressions for the chemical potential 
gives 

Qbb^f — 9bf^b Hri s 
n F = — (10) 

^9bb ~ 9bf 



<±n^ B — yBFH'F 

n B = -rzz -2— (11) 



Att/^b - 9bf^f 
4ng B B -9bf 

The expressions ( TTUj) . (TTTT) provide us with the information on the phase diagram for phase 
separation of the 2D Bose-Fermi mixture. The condition for a finite number of fermions is 



Hf > Qbf^b / 9bb and for a finite number of bosons, it is \lb > QbfI^f / 4vr, where we have 
used 4:7igBB > Qbf^ as required for the stability of the system. To use the above results 
for the BEC side of the BEC/BCS crossover with imbalance, we substitute the bosonic 
chemical potential by hb — ► \i\ + + Ej, = 2/j, + and the fermionic chemical potential 
by Hf —> A*t = A 4 + C- The condition for a finite molecule density is then 

while the condition for a finite fermion density reads 

^ C ~ E b g BF / g BB 

a < : (Id) 

2gBF/gBB - 1 

In between these limits, a mixture of bosonic molecules and fermionic atoms may coexist. 
Thus, we expect corrections beyond mean field to allow for an additional phase to appear in 
the phase diagram of Fig{TJ the imbalanced superfluid, qualitatively indicated by the shaded 
region. Moreover, Eq. ([13]) indicates that the line separating the spin-balanced superfluid 
(molecules only) from the imbalanced superfluid lies below the mean field SF-NP phase 
boundary of Fig{TJ 

In conclusion, the experimental and theoretical investigation of resonance superfluidity 
in Fermi gases has been marked by two recent advances: the creation of a Fermi superfluid 
from a gas with unequal spin populations 0, f| , and the detection of fermionic superfluidity 
in an optical lattice 14| . We have combined both effects and studied the effect of unequal 
spin populations on the superfluidity in an optical lattice. 

We derived an analytical expression for the free energy of the imbalanced 2D Fermi gas, eq. 
([6]). From this result, the BEC/BCS gap equation and saddle-point number equations can 
be solved analytically, and results are provided for fixed imbalance Sn in table (JTJ) However, 
analyzing the free energy as a function of the superfluid gap we found that phase separation 
occurs, and a (possibly imbalanced) normal and (always balanced) superfluid phase separate. 
The phase boundaries are evaluated analytically. The shell structure of the inhomogeneously 
trapped gas is revealed to successively be from the center to the edges a balanced superfluid, 
the spin-imbalanced normal state and the fully polarized normal state. In the deep BEC 
limit a Bose-Fermi mixture offers a more suitable description than the mean-field picture. 
In this Bose-Fermi mixture, an imbalanced superfluid mixture of bosonic molecules and 
fermionic atoms is possible in addition to the phases identified earlier. 
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